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ABSTRACT 

The coronal magnetic field above a particular photospheric region will vanish 
at a certain number of points, called null points. These points can be found 
directly in a potential field extrapolation or their density can be estimated from 
Fourier spectrum of the magnetogram. The spectral estimate, which assumes 
that the extrapolated field is random, homogeneous and has Gaussian statistics, 
is found here to be relatively accurate for quiet Sun magnetograms from SOHO's 
MDI. The majority of null points occur at low altitudes, and their distribution 
is dictated by high wavenumbers in the Fourier spectrum. This portion of the 
spectrum is affected by Poisson noise, and as many as five-sixths of null points 
identified from a direct extrapolation can be attributed to noise. The null dis- 
tribution above 1500 km is found to depend on wavelengths that are reliably 
measured by MDI in either its low-resolution or high-resolution mode. After cor- 
recting the spectrum to remove white noise and compensate for the modulation 
transfer function we find that a potential field extrapolation contains, on average, 
one magnetic null point, with altitude greater than 1.5 Mm, above every 322 Mm 2 
patch of quiet Sun. Analysis of 562 quiet Sun magnetograms spanning the two 
latest solar minima shows that the null point density is relatively constant with 
roughly 10% day-to-day variation. At heights above 1.5 Mm, the null point den- 
sity decreases approximately as the inverse cube of height. The photospheric field 
in the quiet Sun is well approximated as that from discrete elements with mean 
flux = 1.0 x 10 19 Mx distributed randomly with density n = 0.007Mm~ 2 . 

Subject headings: MHD — Sun: corona — Sun: magnetic fields 
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Introduction 



Models of magnetic reconnection, from the very earliest, have assig ned particular signifi- 



cance to points where the magnetic field vanishes: magnetic null points (jDungeylll958l ; [Sweet 



19581 ). In two-dimensional reconnection models, null points (also called neutral points) are 
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Dissipation of magnetic energy at magnetic null points is therefore a potential source of heat 
in the quiet Sun corona. The effectiveness of this dissipation for heating the entire corona 
will depend, in part, in how null points are distributed throughout its field. 

Coronal jets are taken as another observational manifestation of coronal nul l points. Ob- 
serva t ions in EUV and X-r ays show that jets, with their characteristic apex cusp (jShibata et al 
19921 ; IShimojo et q/.lll996f). appear to be more common in coron al holes and the quiet Sun 
than previously believed (ICulhane et a/.ll2007uCirtain et a/.ll2007l ). Theoretical models of jets 
invoke magnetic reconnection occurring at a null point many megameters above the photo 



spher ic surface ( lYokoyama and Shibatal Il996l ; iMoreno-Insertis. Galsgaard and Ugarte-Urra 



20081 ). In order that such a model apply to the large number of jets observed it would be 
necessary that magnetic null points are relatively common at the altitudes of observed cusps. 

In many models, coronal null points occur above a photospheric field consisting of one 



magnetic element completely surrounded by a rin g of opposing polarity (lAntiochosi 11998 



Moreno-Insertis. Galsgaard and Ugarte-Urrall2008l ). The null point will persist even if the 
continuous surrounding ring is broken, but it remains necessary that opposing polarity be 
found on "all sides" of the central element (ISeehaferl 1 19861 ; iBeveridge. Priest and Brown 
20041 ). This suggests that coronal null points occur only under special circumstances which 
are likely to be found infrequently in the actual solar photosphere, and raises once more the 
question of how common null points might be in the actual coronal magnetic field. 

One recent investigation by Regnier, Parnell and Haynes (2008) found 80 magnetic null 
points over a 102 Mm x 116 Mm patch of quiet Sun. A potential field was extrapolated onto 
a rectilinear grid from a magnetogram made by Hinode's NFI instrument. The algorithm of 
Haynes and Parnell (2007) was then applied to the gridded field to locate all points where 
an interpolated field would vanish. This revealed that, at least over that region at that time, 
there was a column of N n = 6.7 x 10~ 3 null points per square megameter, about half above 
a height of z = 1 Mm. Close, Parnell and Priest (2004) extrapolating from high-resolution 
MDI magnetograms of a different day found 1.7 x 10~ 3 Mm -2 , almost four times fewer. This 



- 3- 



discrepancy suggests a potential dependence on instrument and specific date of observation. 
A more general measurement of null point density, applicable to any time or any region, 
would require the repetition of this process. 

Several investigations have estimated the density of magnetic null points in generic quiet 
Sun coronal fields. Schrijver and Title (2002) modeled the photospheric field of the quiet 
Sun as a balanced, random distribution of magnetic point charges. They found that the 
potential field anchored to these charges vanished at one point in the corona, on average, 
for every eleven point charges on the photosphere. They estimated that the density of null 
points diminished exponentially with height. 

Longcope, Brown and Priest (2003, hereafter LBP03) developed a more general approach 
to estimating the density of coronal null points in a potential coronal field. They assumed 
the corona had a random, Gaussian magnetic field, anchored to a statistically homogeneous 
photospheric field with a known power spectrum S(k). The density of null points at a given 
height, Pn(z), scaled with the inverse cube of the characteristic scale length at that height. 
This was expressed in terms of integrals of the power spectrum, S(k), as well as on the 
mean of the vertical field, B z . At lower heights it depended on larger wave numbers (smaller 
scales). A perfectly flat spectrum, such as from ideal white noise or a distribution of point 
magnetic charges, introduces no length scale of its own. In this case the density decreases 
with height as 

0.040 

independent of the strength or other properties of the photospheric field. 

The spectrum from a magnetogram can be used in the spectral estimate of LBP03 to 
compute the distribution of nulls above observed regions. Its assumption of randomness 
and homogeneity are likely to be fairly valid in quiet Sun regions. In this work we quantify 
the degree to which a potential field extrapolation satisfies the other assumption, that the 
field components are Gaussian random variables. We also test the estimate against the 
distribution of null points found directly in an extrapolated field. We find evidence that the 
spectral estimate gives a reasonable approximation of the number of magnetic null points 
above a region of quiet Sun, especially above 1,500 km from the magnetogram surface. 

The spectra of MDI magnetograms reveal several factors which could affect the number 
of null points found either directly or though spectral estimate. There is a level of white noise 
which, while reduced by averaging, is present in any magnetogram. We find that removing 
its contribution from a spectrum reduces the number of null points by as much as 85%. 
The defocusing of the telescope, partly on purpose in the low-resolution mode, steepens the 
spectrum, causing an underestimation of the number of null points. Both of these effects 
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can be removed from the power spectrum, S(k), yielding a spectral estimate of null points 
less affected by instrumental artifacts. 

We perform spectral estimates on 562 MDI magnetograms of quiet Sun during the two 
solar minima it has observed (1996-1998 and 2006-2007). After compensating for the effects 
mentioned above, we find the number of coronal null points above 1.5 Mm to be fairly 
constant: N n = 3.1 x 10~ 3 ± 3.0 x 10 _4 Mm -2 That is, on average, one such null point over 
every 322 Mm 2 patch of quiet Sun. When the Sun is completely free of active regions this 
amounts to about 20, 000 null points above the solar surface. 

The magnetic spectrum of the quiet Sun is found to be remarkably consistent over 
the periods investigated. The typical spectrum is consistent with a potential field from 
magnetic point sources located about 1.4 Mm below the photospheric surface. The variance 
and kurtosis of the magnetograms can be used to infer the size and areal density of these 
sources. Assuming them to have an exponential distribution of fluxes we find the sources 
have mean flux, = 1.0 x 10 19 Mx and areal density n = 0.007 Mm -2 . 

The paper is structured as follows. In the next section we assess the degree to which po- 
tential fields, extrapolated from two particular MDI magnetograms, satisfy the assumptions 
of the LBP03 spectral estimate. We then directly locate all of the null points in extrapo- 
lations from these magnetograms and compare the results to the spectral estimate. In §3 
we apply the spectral estimate to 562 MDI magnetograms spanning solar minimum periods 
in 1996-1998 and 2006-2007. We find that the number of coronal null points stays at a 
relatively constant level, with ~ 10% day-to-day variation. In §4 we show how the generic 
spectrum of the quiet Sun photosphere may be interpreted as a random super-position of 
identical magnetic elements. We infer the properties of the elements and their distribution 
from the Fourier spectrum. Once again, these properties appear relatively constant over 
time. Finally, we discuss what our findings imply about null points in the actual coronal 
field above quiet Sun. 



2. Distribution of Nulls 

2.1. The Gaussianity of the extrapolated field 

The basis for our computation will be magnetograms from regions of quiet Sun, near 
disk center. Figure [1] shows two examples of such magnetograms, each made by averaging 
five one-minute cadence MDI line-of-sight magnetograms. The field of view is 300" x 300" 
centered at disk center. The left and right images are from high-resolution and low-resolution 
data obtained 1.75 hours apart. 
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Fig. 1. — Magnetograms of a quiet sun area from 23 Feb. 1997, made by averaging five 
consecutive MDI magnetograms made at one-minute cadence The grey scale shows the line 
of sight field between ±100 Gauss, (a) High resolution magnetograms and (b) Low resolution 
magnetograms from 1:45 later. 



The spectral estimate of LBP03 yields the spatial density of null points in a coronal 
magnetic field satisfying certain conditions. It assumes the coronal field to be potential 
and to be anchored to a random photospheric field whose statistics are homogeneous and 
Gaussian. These conditions on the photospheric field (the magnetogram) lead to a coronal 
field whose components are statistically homogeneous and Gaussian. From the properties of 
this field LBP03 reports the mean density of null points at any height in the corona. 

Extrapolating upward from the magnetogram of Figuredh to a particular height z results 
in a new field, B z (x,y,z), whose histogram we find to be nearly Gaussian. Figure [2] shows 
histograms from several heights as solid lines, and the Gaussians best fitting their central 
portions as dashed curves. The cumulative distribution from a range of values was fitted, 
and the range was iteratively adjusted until it contained exactly two standard deviations to 
either side. The standard deviation of this fit, ai 9 \ is indicated by triangles and plotted on 
axes to the left. The standard deviation of the entire field, ai m \ plotted with a dashed line 
on the same axis is systematically greater than that of the core. 

It can be seen in Figure [2], that the field B z has an approximately Gaussian distribution 
even well above the photospheric level (z = 0). At the bottom, the Gaussian component has 
a standard deviation, ai g ^ ~ 11 .3 Gauss, comparable to the intrinsic noise in a five- minute 



average of MDI magnetograms (ILiu and Nortonll200ll ); it decreases with height as the width 



of the Gaussian narrows. The actual distribution departs from the Gaussian fit principally 
by an excess (tail) mostly outside ±ai g \ Histograms made from the horizontal components 
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Fig. 2. — Histograms of the field B z computed by extrapolating the photospheric magne- 
tograms to five different heights. The height, in Mm, of the slice corresponds to the level of 
the bottom of the histogram on the far left axis. A Gaussian fit to the central region is shown 
as a dashed line, with triangles at ±ai 9 \ These values are also plotted against a horizontal 
axis on the left. A solid curve shows the full functional form of ai 9 \z). A dashed line shows 
(ri m \z) on the same axis (squares occur at heights coincident with the histograms). 
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of the extrapolated fields at varying heights look very similar to those of the vertical field 
shown in Figure [2j 

The Gaussian fit to the central core of each distribution can be used to quantify the 
degree to which the actual field is a Gaussian variable. The amplitude of the Gaussian 
component ranges from 0.97 at z = to 0.81 higher up (see Figure [3b). That is to say the 
field is between 81% and 97% Gaussian. 




Fig. 3. — (a) Standard deviations from all three components of the potential extrapolation of 
the magnetogram Figure [Ik. The solid and dashed curves show and respectively for 
B x (crosses), B y (triangles) and B z (squares), plotted versus height. The symbols occur at 
the same heights as the histograms from Figure [2j Dotted lines show the standard deviation 
for vertical (upper) and horizontal (lower) field components extrapolated from white noise, 
(b) The fraction of each field which is fitted by the Gaussian within the central Aa^. 

While a significant component of the magnetic field has a Gaussian distribution, it 
is not spatial white noise. Longcope, Brown and Priest (LBP03) show that a potential 
field extrapolated from a magnetogram with isotropic power spectral density (psd) S(k) has 
vertical field with variance 

POO 

a 2 z (z) = 2tt / S(k)e- 2kz kdk = 2a 2 x (z) = 2a 2 (z) . (2) 
Jo 
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Spatial white noise is uncorrelated from pixel to pixel and therefore has a flat psd 



S(k) 




e(k c -k) 



(3) 



out to the resolution limit k c = 2 v / 7r/ Ax for pixel size Ax (0 is the Heaviside function). 
Using this in (T5]) gives the variance due to spatial white noise alone 



This functional form, plotted as a dotted line for each component of the field in Figure [31 falls 
off much faster than either or . We therefore conclude that the Gaussian component 
consists of contributions with spatial correlations longer than a single pixel: it is not white 
noise. 

It is common practice to assess the degree to which a variable departs from Gaussian 
by its kurtosis statistic, K. The kurtosis is the fourth moment divided by the square of 
the variance, less three; when the variable is perfectly Gaussian K — 0. Excesses, relative 
to a Gaussian, in the tails of the distribution will elevate the fourth moment and lead to 
a positive kurtosis. Figure H] shows that even when the excess tails account for 3% of the 
distribution, they lead to a large kurtosis. 



Fig. 4. — The kurtosis of the fields for B x (crosses), B y (triangles) and B z (squares), plotted 
versus height. Solid and dashed lines are for the high resolution and low resolution magne- 
tograms respectively. The dotted lines show the results for point magnetic sources, taken 
from Equation (J6j). 

For comparison, a distribution of n point sources per unit surface area has a flat power 
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spectrum with no cut-off (k c — > oo) whose variance is (LBP03) 

«JM = ^ , (5) 



where (</2 2 ) is the variance in the source fluxes. The kurtosis of this same photospheric field 
is 

= Ai^A^ = 4^ = 4K V . (6) 

5tt (^) 2 nz 2 x v K J 

Both of the expressions above diverge at the photosphere (z — > 0) as a consequence of the 
very strong fields in the neighborhood of a point source. This divergence is one element 
leading to the failure of the spectral estimate at low heights, z < rr x l 2 (LBP03). 

It is evident from Figure H] that the kurtosis of the actual field is generally large but 
remains finite even as z — > 0. The large magnitude arises from the tails in the pdf corre- 
sponding to those pixels of the magnetogram in magnetic elements, where \B Z \ > 5a 2 (0) ~ 50 
Gauss. These are the excess tails of the distribution function, but they are also the structures 
which generate field above a few megameters. 



2.2. Assessing the accuracy of the spectral estimate 

A spectral estimate of the magnetic null point density is found from the isotropic power 
spectrum, S(k), of the photospheric magnetogram. The spectrum determines the variance, 
as given by Equation (j2J), and a related quantity 

q 2 (z) = S(k)e- 2kz k 3 dk . (7) 

°z\ Z ) JO 

Provided the field has Gaussian statistics the density of null points in it is found to be 

p N (z) = G{B z /a z )q\z) , (8) 

where G(() is a dimensionless function defined in LBP03. The function depends on the mean 
vertical field strength B z and also weakly on another moment of the spectrum. Provided the 
polarities are sufficiently balanced that B z a z over the relevant range of heights, we find 
G~2x 1CT 2 for most spectra of quiet Sun magnetograms. 

We assess the accuracy of the approximation offered by the spectral estimate by direct 
comparison. Both magnetograms from Figured] are extrapolated onto a rectilinear grid with 
periodic lateral boundaries. Between the grid points the magnetic field is defined using 
tri-linear interpolation. A magnetic null point is the point where all three components of 
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this interpolated field vanish together. These null points are located using the algorithm 
of Haynes and Parnell (2007). The number of null points found above a given height z is 
divided by the photospheric area to yield a null-point column density, N n (z). The right and 
left axes of Figure [5b shows the number and column densities found this way, respectively. 
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Fig. 5. — The density (a) and column density (b) of null points in potential fields extrapolated 
from the MDI magnetograms in Figure HJ Dark and light curves correspond to the high 
resolution and low resolution data respectively. The solid curves are computed by finding 
the nulls on a gridded extrapolation; the dashed curves are the spectral estimates. Symbols 
appear at heights of one, two, three and four horizontal pixel sizes (z = Ax, 2Ax, 3Ax and 
4Ax). For panel (b) the left axis reports the column density in nulls per Mm 2 above a given 
height, while the right axis gives the total number over the entire magnetogram. 

There are more null points in the high resolution case (dark curves) than in the low 
resolution (6135 vs. 664). The majority of the excess null points lie in a layer about one Mm 
thick. Over heights above one Mm the predicted null columns (dashed curves on Figure [5b) 
from the high resolution and low resolution cases are quite similar. The greatest disagreement 
between the actual null column, N n (z), and the spectral estimate occurs within this same 
one Mm bottom layer. Above this height the two estimates are within a factor of two of 
reality for both resolutions. Since the bottom layer of the high resolution field has the largest 
kurtosis it is perhaps unsurprising that the greatest disagreement with the spectral estimate 
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occurs there. Perhaps more surprising is that kurtosis as large as 40 (the maximum of the 
low resolution data) leads to disagreement by not more than a factor of two. 



3. The quiet Sun null column in general 

3.1. Comparing nulls from different magnetograms 

The extent to which, for levels above one Mm, the spectral estimates of high resolution 
and low resolution agree can be understood from a comparison of their power spectra S(k), 
as shown in Figure El A fast Fourier transform (FFT) of the L x x L y magnetogram returns 
an array of complex coefficients, -B m> „, for the expansion 

B z (x,y) = B m>n exp[ % mAk x x + i nAk y y ] , (9) 

m,n 

where Ak x = 2tt/L x and Ak y = 2ir/L y . The quiet-Sun magnetograms have very few fea- 
tures at their edges which might suffer discontin uity under period ic repetition so we find it 



unnecessary to employ any windowing function (IPress et al. 



19861 ). 



The two-dimensional power spectrum at wave number k = mA^x + nAk y y is 



\B 



12 



s < k > - Hkf, ■ < 10 > 

This function is defined only on points of the Fourier-space grid, but approximates a contin- 
uous function. The denominator in Equation (fTO]) eliminates from S*(k) any dependance on 
pixel size or field of view. 

The isotropized power spectrum, S(k), is computed by averaging all values of S(k) 
within an annulus of Fourier-space. We use annuli of radii extending to a maximum wave 
number k c = 2y / 7r/Ax, chosen so that the total Fourier-space area, 7fk 2 , matches that 
occupied by the actual magnetogram. The cut-offs for high and low resolution MDI magne- 
tograms are therefore 8.2 rad/Mm and 2.5 rad/Mm respectively since these have pixel sizes 
of Ax = 0.43 Mm and Ax = 1.42 Mm in this particular magnetogram. 

The spectrum in Equation fflOl) is normalized so that, 

(B 2 z ) = ^|5 m>n | 2 ~ J S(k)d 2 k = 2n C S(k) kdk , 



rn.n 



where the angled brackets are a spatial average. The variance is found from the same 
expression but excluding from the sum (m,n) = (0,0). The variance at height z above the 
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Fig. 6. — The isotropized spectra for 300" x 300" magnetograms of quiet Sun centered at disk 
center on 23 Feb. 1997. Curves of different colors are S(k) for 30-second MDI exposures at 
high (green) and low (red) resolution; 5-minute averages at high (magenta) and low (black) 
resolution; and Kitt Peak magnetogram (blue). The inset is an expansion of the region at 
low wavenumbers. 
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magnetogram plane, Equation 
extrapolation. 



includes an exponential factor from the potential field 



Each spectra from an MDI magnetogram in Figure M has a range of white noise over 
its highest wave numbers. We estimate its level, S n , from the median value of S(k) over the 
highest 15 spectral bins. The noise level is reduced by exactly 1/5 in the a yerages strongly 
suggesting an origin in Poiss on counting noise, uncorrelated between frames (]Liu and Norton 
200ll ; beForest et~al\ I2007T ). The fiat noise floors contribute a 2 (0) = irk c S% = (8.8 G) 2 
and (6.3 G) 2 to the low and high-resolution variances respectively. We noise-correct each 
spectrum by subtracting from the entire curve the flat noise floor, S n , determined from the 
top 15 spectral bins. 

Each of the noise-corrected spectra appears to hav e been steepened by t he modulation 
transfer function (MTF), M(k), of its imaging system (IBorn and Wolfl Il980l ) . and must be 
corrected to S(k)/M 2 (k). The most conservative correction to the high-resolution spectrum 
is to assume it to be diffraction-limited and take M(k) = Md(k), given in an appendix. MDI 
is intentionally defocussed in its low-resolution mode, leading to an MTF we approximate 
as M 2 (k) = M%(k)/(1 + £q/c 4 ) where £q ~ 1.6" at this time (see appendix). Performing 
these corrections restores all four MDI spectra to similar forms, over low wave-numbers, 
as shown in Figure [7J The Kitt Peak spectrum (blue) remains below all others since it 
has not been divided by its MTF, part of which pres umably arises from atmospheric seeing 
(IHufnagel and Stanley! 1 19641 ; lAbramenko et all l200ll ) . 



There is, over the range k > 1.5rad/Mm, an excess in high-resolution power, relative 
to low-resolution; this excess is further amplified by the MTF correction. The excess might 
be partly an artifact f rom the granular mot ions during the brief intervals between different 
filter-gram exposures (IDeForest et al\ 120071 ). Its origin, whatever that may be, is evidently 
correlated several-fold longer than five minutes, since averaging diminishes it by less than a 
factor of five. We next show that this power does not significantly affect the null densities 
at heights where potential field extrapolation is warranted. 

The spectra from Figure O are used in the Laplace transform integrals, ([2]) and (j7J), 
from which are computed the null density Pn( z ), according to Equation (jSJ). This density 
is integrated downward from some upper height yielding the null column density, N n (z), 
plotted in Figure [HJ The exponential factor, e~ 2hz , in each of the spectral integrals means 
that the null density at a given height z depends principally on the spectrum over wave 
numbers k < 1/z. Thus heights above z = 1 Mm depend mostly on the spectrum below (to 
the left of) k = lrad/Mm. As discussed above, and evident in Figure UJ all MDI spectra 
are fairly similar over this range. It is therefore natural that all five solid curves in Figured 
follow similar paths above one Mm. 
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Fig. 7. — The same spectra from Figure [6] after removal of white noise and correction for 
the MTF. No corrections are performed on the Kitt Peak spectra (blue) which is retained 
for reference. Dashed curves show the high resolution spectra before it was corrected for the 
MTF, but after white noise was removed. 
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Fig. 8. — Spectral estimates of the null column densities, N n (z), computed from each of the 
spectra in Figures [6] and [7J using the same color code. Dashed curves are from the original 
spectra (Figure [7j); solid curves are from spectra from which white noise has been subtracted 
and which have been divided by the MTF (Figure [7]). Dotted lines show the curves after 
noise subtraction without MTF correction. 
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Noise contributes significant power to the highest wavenumbers, creating a population 
of artificial null points in an extrapolation. The null column densities after subtraction of 
this noise, the dotted curves in Figure El lie below the corresponding uncorrected (dashed) 
curves. The fraction of nulls attributable to noise alone is lower in five-minute averages than 
one-minute magnetograms, since the former have lower noise. The effect is most pronounced 
in one-minute low-resolution magnetograms, for which 85% of all nulls are due to noise. 

The small-scale structure revealed in high resolution magnetograms translates into com- 
plex fields with high null densities at the lowest level of the potential magnetic field. While 
potential fields are frequently used to approximate the quiet Sun corona, they are of ques- 
tionable validity in the lowest layers, where the plasma density and pressure are liable to 
play their most significant role. Above a height of 1,500 km (1.5 Mm) the potential field ap- 
proximation may be on better footing. It appears from Figure [8] that most of the estimates 
based on potential field extrapolation, are in reasonable agreement above this particular 
height (at z — 1 Mm, for example, the high and low resolution estimates appear to disagree 
to a greater extent). For these reasons we hereafter focus on estimating the total column of 
null points above z = 1.5 Mm, denoted iV„(1.5Mm). 



3.2. The generic behavior 

The foregoing sections have considered different magnetograms of the same region of 
quiet Sun (23 Feb. 1997). In order to generalize our conclusions we perform similar analysis 
on magnetograms spanning solar minimum from 2006 and 2007. The same central 300" x 300" 
sub-image is extracted from one full-disk (low resolution) daily MDI magnetogram every 
second day. It is analyzed provided there is no obvious active region field, or preponderance 
of one polarity, and provided there is less than 8 x 10 19 Mx in pixels above 200 G. This 
value was selected after some trial and error to eliminate plage regions whose presence would 
violate our assumption of homogeneity. Separate analyses are performed for one-minute 
daily magnetograms and for those composed of 5-minute averages. 

The isotropized spectra from all quiet Sun areas turn out to be very similar to one 
another, closely resembling the example from the previous section. Figure [9] shows a median 
spectrum: the the median value, at each wave-number, from over 110 quiet Sun spectra of 
each low-resolution kind. The upper and lower quartile points, bracketing half of all values, 
are plotted as lighter curves. This clearly shows that the most significant variation occurs 
at wavenumbers below 0.2 rad/Mm (i.e. wavelengths longer than 30 Mm). 

Removing the white noise floor by first subtracting it and then extrapolating the noisier 



0.0 0.5 1.0 1.5 2.0 2.5 

k [rad/Mm] 

Fig. 9. — Median spectra from magnetograms in 2006 and 2007. (a) The median spectrum 
(dark) and the upper and lower quartiles (light curves) from the five-minute averages, (b) 
The median from 118 one-minute magnetograms (upper curve) and 114 5-minute averages 
(lower curve). The lighter curves show the results of subtracting white noise (lower pair) 
and then dividing by the MTF (upper pair). 
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portions yields the lower light spectra in Figure [9b. This is then divided by the empirically 
determined MTF, M 2 (k) = Mj(k)/(l+£fo 4 ), where £ = 1.04" as described in the appendix. 
The results are the two upper light curves. 

The column density found from these various versions of the median spectra are shown 
in Figure [10J The solid curves of each shade are the result of removing noise and correcting 
for the MTF. The broken curves above these are from the raw spectra, while the dashed 
curves below are after noise has been removed but no MTF correction performed. Symbols 
mark the values of these at z — 1.5 Mm: 3.0 x 10~ 3 Mm -2 and 2.9 x 10~ 3 Mm -2 respectively. 




2 4 6 

height z [ Mm ] 



Fig. 10. — The column densities from each of the median spectra in Figure [9b. Dark and 
light curves are for one-minute magnetograms and five-minute averages respectively. Solid 
curves are from the entire spectra and dashed curves are from the noise-corrected spectra. 
A dotted curve is an empirical fit. 



The median spectra are clearly very close to pure exponentials, S(k) ~ e~ 2kd . This 
functional form suggests the empirical model for the null column density of one-minute 
magnetograms (dotted line) 

AM*) = ^ , (ID 

where the numerator and d = 1.6 Mm are found from fitting. 
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The null column above a given height does vary between magnetograms, as shown 
in Figure [TTJ There appear to be random fluctuations superimposed on a much slower 
variation. We estimated the dispersion about the trend, 2.7 x 10~ 4 Mm -2 for the one-minute 
magnetograms and 2.5 x 10~ 4 Mm -2 for five- minute averages, by subtracting a 90-day running 
mean. Considering the entire data set of one-minute magnetograms, including its slow trend, 
we find the column of null points above z = 1.5 Mm to be 

iV n (1.5Mm) = (3.1 ±0.3) x 10~ 3 Mm~ 2 . 

(The errors quoted above represent rms variation; the error in the mean is far smaller.) The 
value from five-minute mae;netoe;rams is slightly lower: (2.9 ± 0.3) x 10" 3 . We believe this 
value is lower partly because of motional blurring during the averaging interval, so we adopt 
the one-minute-magnetogram value as the more accurate of the two. 
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Fig. 11. — The column of null points above z = 1.5 Mm found from quiet sun magnetograms 
spanning 2 years of solar minimum. Triangles and crosses denote one-minute and five-minute 
low-resolution magnetograms respectively; squares are from high-resolution magnetograms. 
Error bars on the right show one standard deviation of scatter about a 90-day running mean. 
Solid lines along the bottom show the reported focus position of the MDI telescope. 

We believe the slow variation in iV n (1.5) is due to gradual changes in the MTF of 
the imaging system. The reported focus position of the instrument is shown as a solid 
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curve long the bottom. The focus is periodically changed to compensate for changes in the 
instrument itself. Our estimate of the MTF was performed through cross-comparison to the 
high resolution images, plotted as squares, all while the instrument was in focus position 
5. This is the MTF by which we are correcting the spectrum, so we believe that the value 
in the middle of this interval is probably the most accurate. The MTF correction raises 
the low resolution columns by about 40%, thereby making them more consistent with high 
resolution values. Omitting that step gives null columns which could be considered lower 
bounds on the actual values: medians are iV n (1.5) = 2.2 x 10- 3 Mm" 2 for one-minute and 
2.1 x 10~ 3 Mm~ 2 for five-minute. 

A comparable plot from the previous solar minimum, Figure [T2l shows more pronounced 
evolution on long time scales. This data comes from the very beginning of the SOHO mission, 
when the instrument was probably changing most rapidly. There is a very large change in 
Nov. 1997 coincident with the refocussing of the telescope. The later refocussing events 
in Figure [11] produce much less obvious changes, either because they are more frequent or 
because there is less dramatic change between them. Due to the larger, and more rapid 
change we did not attempt an estimate of the MTF. Instead we used the empirical form 
from the later phase, with £ = 1.6" found from fitting the 23 Feb. 1997 data described 
above. This restores the typical values to the same range as they were in the 2006/2007 data 
set. 



Some insight into the origins of the quiet sun magnetic field may be gleaned from 
its spectrum. Consider a model photospheric field composed of distinct elements, each 
characterized by the same normalized shape function, /(x). Elements at locations, Xj = 
(xj,yj), with fluxes, <pj, will create a vertical field 



in the photospheric plane. For simplicity let us assume that the fluxes are random, have 
zero mean, and are uncorrelated: {4>i4>j) = {4> 2 )^ij- We also assume that the fluxes are 
independent of their locations x, which are distributed randomly with areal density n. Under 
these assumptions the power spectrum of the photospheric field is 



4. Modeling the quiet Sun spectrum 



4.1. A synthetic magnetogram 
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Fig. 12. — Null columns during the solar minimum from 1996-1998. Plot format is the same 
as in Figure [TT1 



where /(k) is the Fourier transform of the shape function. 



A purely exponential spectrum can be interpreted as being composed of individual flux 
elements each with an exponential Fourier transform such as 



This function also happens to be the vertical potential field due to a unit point charge located 
at z = — d; it is the potential field Green's function. 

The observed spectrum, before MTF-correction, is very close to an exponential, S(k) ~ 
e~ 2kd . The exponent is thus directly related to a "depth" at which point sources are sub- 
merged. On the other hand, the spectrum after correction, S(k) ~ e~ 2kd /M 2 (k), has an 
inverse Fourier transform which is not a simple function, and does not admit such a sugges- 
tive interpretation. Figure [13] shows the functions f(r) for both the uncorrected (dashed) 
and MTF-corrected (solid) spectra from 23 Feb. 1997. The integrated field (bottom) shows 
that for each case half of all flux is confined to a similar tq ~ 4 Mm circle. 

The variance of a photospheric field composed of random elements is 



where r* characterizes the radius of the shape function, f(r), through the integral of its 
square. The parameter is defined to match the root-mean-squared radius of a Gaussian 
shape. The uncorrected shape function, Equation (1T41) . yields, r* = \f2d, whose introduction 
into Equation ffl5l) returns the variance of the vertical field at a height d above a distribution 
of point charges, i.e. Equation (j3J). The MTF-corrected shape function is notably narrower 
(solid line in Figure [TBI, under the anticipation that it is broadened by the point spread 
function to produce the broader observed shape (dashed) with its larger radius r*. 

The kurtosis of the composite field given by expression (Fl2|) is 



where the dimensionless coefficient x depends only on the shape of the individual elements, 
f(r). It is unity for a perfectly flat function and greater as the shape becomes more peaked. 
For charges submerged at a depth d, the value x = 32/5 can be inferred by comparison 
between expressions flBD and (fTE|) and use of the relation r* = \/2d = \/2z. 

The dimensionless ratio, £ = (0 4 )/(0 2 ) 2 , depends only on the distribution of element 
amplitudes, and not on their size or shape. Flux amplitudes distributed uniformly or expo- 
nentially will be characterized by £ = 9/5 or £ = 6 respectively. Parnell (2002) found that 
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Fig. 13. — The structure of a magnetic element inferred from the power spectrum of a 
magnetogram. The top curves show the shape functions, f(r), from Equation (114j) . for 
an uncorrected exponential spectrum with d = 2.2 Mm (dashed) and for the shape after 
correction for the low-resolution MTF with £ = 1.6" (solid). Symbols on each curve mark 
r*. The vertical scale is the same for each curve and makes 2ir J f(r)rdr = 1. The bottom 
panel shows the total fraction of flux with a given radius, and symbols mark the radius of 
half flux. 
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fluxes of each polarity were distributed with a hybrid between an exponential and power-law 
known as a Weibull distribution. The particular Weibull distributions reported are char- 
acterized by values of £ between 22 and 30. These are larger than for purely exponential 
distribution since broader wings lead to more common occurrence of large fluxes. 

The variance of Gaussian white noise, er 2 , adds to expression (Ti~5]) to give the total 
variance in a magnetogram, cTq = a 2 + a 2 . The Gaussian noise has zero kurtosis, so its 
addition will decrease the kurtosis of the composite, Equation (flBT) . by the factor cr^/a^. 
The areal density of flux elements follows directly from the magnetogram's kurtosis, K , and 
properties of the elements 



Having found the density, expression (Tl5|) can be used to find the root-mean-square (rms) 
amplitude of the flux, (0 2 ) 1 ^ 2 . 

Using the above reasoning we may model a quiet sun magnetogram as a combination 
of Gaussian white noise and a random distribution of flux elements. The magnetogram in 
Figure [lb, for example, has a power spectrum which, after removal of white noise, fits an 
exponential S(k) ~ e~ 2kd , with d ~ 2.2 Mm. The magnetogram has a standard deviation 
o"o — 15.8 G, only a n = 8.8 G of which comes from white noise. Assuming an exponential 
distribution of flux amplitudes (£ = 6) the kurtosis, K = 38, implies a flux element density, 
n = 0.0039 Mm -2 , from Equation (TlTl) . Equation ([15]) then gives the rms flux, (0 2 )^ 2 = 
2.3 x 10 19 Mx. The exponential flux distribution and zero mean value can be used to convert 
the rms to = 1.6 x 10 19 Mx. 

Figure [T4"a shows a magnetogram synthesized from the distribution of sources described 
above. Its parameters were chosen to resemble the five-minute average low resolution magne- 
togram of [lb, and the two have indistinguishable power spectra (Figure [Mb) . The standard 
deviation and kurtosis of the synthetic magnetogram, a = 16.1 G and K = 39, are well 
matched to the original. 

Correcting the spectrum for an MTF with l§ = 1.6" leads to a more narrowly peaked 
element, /(r), whose radius is r* = 2.1 Mm and which has x = 13.2 (compared to r* = 
\pld = 3.1 Mm and \ = 6-4 without correction; see Figure [TBI . These more peaked elements 
must be distributed with greater areal density and smaller mean flux in order to match the 
observed kurtosis and standard deviation: n = 0.016 Mm -2 and = 5.6 x 10 18 Mx. The 
instrument's point spread function blurs each element reducing the total apparent variance. 
The MTF correction thereafter amplifies the observed variance leading us to the inference 
of a higher density of fundamental elements. In this case by a factor of four. 
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Fig. 14. — A magnetogram synthesized from a distribution of point sources at d = 2.2 Mm 
below the photosphere, (a) The magnetogram with grey scale over \B Z \ < 100 G, exactly as 
for[Tb. (b) The power spectra of this magnetogram (dark) and of Figure [Tb (light). 



The same analysis can be performed on all of the spectra analyzed in the previous 
section. The MTF correction implies fundamental elements with a median radius r* = 1.1 
Mm, too small to observed directly in low resolution. As a consequence the median element 



density is quite large, 



n 



0.046 Mm , and median flux quite small, 



0.22 x 10 



19 



Mx. These val ues are outside the range found by direct identification o f flux elements in 
magnetograms (jSchrijver et all 119971 ; iParnelll |2002| ; iDeForest et all 120071 ) since the correc- 
tion for point spread function is not generally performed. Indeed, our final results depend 
quite sensitively on this correction, since unresolved elements contribute only slightly to the 
statistics of the observed magnetic field. 

We repeat our analysis without the MTF correction to obtain an estimate which is at 
once more conservative and more consistent with previous investigation. Figure [15] shows 
the depth (bottom panel), areal density, n (middle panel) and flux, (\(f>\) (top panel) from 
low resolution magnetograms over the 2006-2007 solar minimum. These exhibit the same 
slow variation as the null density, presumably for the same reason: gradual variation of 
focus. The MTF-corrected version (not shown) has still more pronounced variation due to 
its amplification of the high-wavenumber portion of each spectrum. The median depth is 
d = 1.43 Mm, and without the MTF correction the median values for areal density and mean 
flux are n = 0.007Mm" 2 and (\<p\) = 1.0 x 10 19 Mx respectively. These are the value for 
which the observations provide direct evidence and are consistent with studies w hich directly 
identify flux elements ( jSchrijver et a/.lll997l ; |Parnelll 120021 ; IDeForest et al\ 120071 ). 
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Fig. 15. — Parameters for the model magnetograms derived from spectra of one-minute 
(triangles) and five-minute (crosses) full-disk magnetograms. (top) The mean unsigned flux, 
in units of 10 19 Mx. (middle) The density of sources in units of Mm' 1 , (bottom) The 
depth d at which the sources are placed. 
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4.2. Using the model in a noise-free spectral estimate 

Submerged point sources alone, without the addition of white noise, have the same 
purely exponential power spectrum as the noise-corrected magnetograms. It was found in 
LBP03, however, that the spectral estimate grossly overestimates the column of null points 
above a distribution of submerged point sources. This was attributed to a severe departure 
from gaussianity within a layer extending Az ~ n -1 / 2 from the source plane (z < 15 Mm in 
our case) . One symptom of this departure is the divergence of the kurtosis given by Equation 
(J6]). The known overestimation for point sources raises concerns about the accuracy of 
the spectral estimate from the noise-corrected magnetogram spectrum, which has the same 
spectrum. 

To further explore this comparison we reduced the noise in the low-resolution 5-minute 
averaged magnetogram, Figure [lb, using a Fourier filter. Amplitudes of all high-spatial 
frequency Fourier modes (|k| > 0.7rad/arcsec) were rescaled to fit a pure exponential, S(k) oc 
e -2kd w ith ft = 2.2, while retaining their original phases. The result was a magnetogram (see 
Figure [TBI) which is visibly cleaner, albeit blurrier than the original (Figure [lb). The central 
core of its distribution function fits a Gaussian of width ai 9 ^ = 7.0 G. The original, 5-minute 
averaged magnetogram had ai 9) = 10.6 G of which k c ^S^ = 8.8 G was found to be white 
noise, leaving a z = 5.9 G to be be genuine. Thus it would seem our Fourier filtering removes 
a substantial amount of the noise. 




Fig. 16. — The result of Fourier filtering the low-resolution 5-minute average magnetogram 
Figure [Tb. (left) The magnetogram plotted using the same grey scale as the original, (right) 
The null densities, Pn{z), for the Fourier-filtered magnetogram (diamonds) compared to 
the spectral estimate (dashed). The null densities from the synthetic magnetograms created 
with low and high densities of submerged point sources are plotted with triangles and squares 
respectively. 
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The Fourier-filtered magnetogram has, by construction, a psd similar to a magnetogram 
synthesized from submerged point sources (like Figure [H^ without added noise). Such a 
magnetogram was synthesized, extrapolated onto a grid and all its null points located using 
the algorithm of Haynes and Parnell (2007). The result was a column density iV n (1.5 Mm) = 
0.3 x 10 -3 Mm -2 of null points. Repeating the procedure for the Fourier filtered magnetogram 
yielded a column, N n (1.5 Mm) = 3.6 x 10~ 3 Mm -2 , more than an order of magnitude greater. 
The spectral estimates for both cases (which have similar spectra) is iV n (1.5 Mm) = 2.0 x 
10~ 3 Mm~ 2 ; this is smaller than the results from the previous section, which also include 
MCT-correction. 

The spectral estimate exceeds the column above point sources by about a factor of 
seven, in agreement with the finding of LBP03 (they find an actual column N n (—d) = 
0.074n = 0.29 x 10~ 3 Mm~ 2 ). On the other hand, it underestimates the column above the 
filtered magnetogram by 44%, in spite of its very large kurtosis: K z = 49.0 (larger than the 
magnetogram before filtering). 

To explore the significance of kurtosis, a second synthetic magnetogram was constructed 
with a greater areal density (n = 10 _2 Mm -2 ) of point sources at the same depth (d = 2.2) 
but with correspondingly smaller flux ({\<p\) = 0.69 x 10 19 Mx) to endow it with the same 
spectrum as the other two. The kurtosis of the resulting magnetogram (K = 29) is lower 
than either of the others. A gridded extrapolation of this field was found to have a column 
iV n (1.5 Mm) = 1.0 x 10 _3 Mm~ 2 . Moreover, the actual density, Pn( z ) falls well below the 
spectral estimate below z = n -1 ' 2 — d = 7.8 Mm, as shown in the right panel of Figure [TBI In 
contrast, the null density of the Fourier-filtered magnetogram tracks the spectral estimate, 
but is slightly above it, over the entire range of heights. The height at which the null density 
of the first synthetic magnetogram becomes accurate, is z = n~ l l 2 — d = 13.8 Mm, consistent 
with Figure [T5b. 

These experiments suggest that the accuracy of the spectral estimate in predicting the 
actual null point column, is unaffected by the absence of white noise. For both the filtered 
and unfiltered 5-minute magnetogram it underestimates the column by a similar factor. This 
performance is far better than it would be if the photospheric field consisted of submerged 
point sources, as its spectrum suggests. The discrepancy cannot be explained by kurtosis of 
the fields, since the spectral estimate is still too large for a point-source-derived field with 
smaller kurtosis. It seems that the spectral estimate owes its success in the quiet Sun to a 
field intrinsically more complex than one from point sources. 
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5. Discussion 

The foregoing has shown that the Fourier spectra of quiet Sun magnetograms have an 
invariant basic form, at least over moderately high wave numbers (0.2rad/arcsec < k < 
1.5 rad/arcsec). This spectral range plays the greatest role in dictating the distribution 
of magnetic null points in the corona above. As a consequence we find that the density 
of magnetic null points in a potential field extrapolation is about JV„(1.5Mm) = 3.1 x 
10~ 3 Mm~ 2 over any patch of quiet Sun, varying by roughly 10% from day to day. This 
value is just more than twice the number of nulls found by Regnier, Parnell and Haynes 
(2008) at heights greater than 1.5 Mm, using Hinode data of much higher resolution. The 
discrepancy could possibly arise from the much greater flux imbalance over that particular 
magnetogram. The colu mn we find corresponds to 2 nulls ab ove every 5 hexagonal, 14-Mm- 



diameter supergranules (IHagenaar. Schrijver and Title! 119971 ) . 



The potential magnetic field extrapolated from a photospheric magnetogram tends to 
become spatially smoother at increasing heights. This is a general property of harmonic func- 
tions and their gradients. One consequence is that the density of null points in such a field 
decreases with height and most of the null points are found close to the surface. Of the null 
points above z = 1.5 Mm roughly 75% are found within the next 3.0 Mm. We find the spec- 
tral estimate of the null point density in the generic quiet Sun to be p n (z) ~ 0.04/ [z + d)~ 3 , 
where d ~ 1.4 Mm. This simple functional form follows from the approximately exponential 
nature of the photospheric spectrum. 

The density at very low heights, say z < 1 Mm, is dictated by very small horizon- 
tal structure in the photospheric field. This structure is the most challenging to observe, 
and its observation is subject to the most instrumental distortion. Moreover, this lowest 
atmospheric layer, the chromosphere, has the greatest plasma pressure and mass density 
and is least likely to satisfy the conditions which would render a magnetic field potential, 
or even force free. For these reasons we choose to focus our investigation on null points 
above z = 1.5 Mm. This nominal height has nothing to do with the empirical depth, d, 
characterizing the exponential slope of the photospheric spectra. Rather, it arose from our 
experience with different instrumentation (MDI low-resolution versus high-resolution) and 
different methods of accommodating instrumental effects such as the MTF. We found that 
the null point density, Pn{ z ), varied less above z — 1.5 Mm due to these instrumental factors, 
so we deemed it to be more accurately measured there. 

The density and distribution of null points in the corona is predicted somewhat reliably, 
although not perfectly, by spectral estimate of Longcope, Brown and Priest (2003, LBP03). 
Magnetic null points can be found directly, and thus with more accuracy, in a potential 
field extrapolated from a photospheric magnetogram. We used this technique to test the 
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spectral estimate, and found it predicted the shape of the density, Pn(z), reasonably well. 
The spectral estimate of the total column was, however, too low by a factor of two to three 
above one Mm. If the spectral estimate for every magnetogram is low by the same factor, 
then the actual null density in the quiet Sun would be greater than the value we report. 

This discrepancy is most likely due to a departure of the actual field from the statistical 
properties assumed by the spectral estimate (homogeneity and Gaussianity) . We did not 
find a satisfactory method of quantifying these factors. One characteristic in particular, the 
kurtosis of the distribution, did not accurately predict when the assumptions were adequately 
satisfied. We did find, however, that the conditions appear to be better satisfied, and the 
predictions better, for actual magnetograms than for a random distribution of point magnetic 
sources. In the latter case the spectral estimate exceeds the number of null points by at least 
an order of magnitude. 

Although it is more accurate, direct enumeration of magnetic null points suffers from 
several disadvantages relative to the spectral estimate. It is more computationally expensive, 
so performing the search on over 500 different grids would be prohibitive. More significantly, 
instrumental artifacts such as noise and MTF are much more difficult to remove from the 
extrapolation than to remove from the spectrum. We found that white noise present in MDI 
magnetograms can contribute artificial null points outnumbering the real ones by as many 
as five-to-one. The ease with which such factors can be removed the Fourier spectra makes 
the spectral estimate a superior one, even given its inaccuracy. 

We found that the spatial spectrum of the quiet Sun magnetic field is well approximated 
as an exponential, at least for k > 0.2 rad/Mm. This spectrum is consistent with the potential 
field from point magnetic sources submerged to a depth of d ~ 1.4 Mm. These manifest 
themselves on the photosphere as circular concentrations of flux, approximately r* = 2.5 Mm 
in radius. We are able to infer the size and areal density of these structures, (\(f>\) = 1.0 x 10 19 
Mx and n = 0.007 Mm -2 , solely from the statistics of the photospheric field. The resolution 
of such structures into finer elements, clearly seen in high-resolution observations, will affect 
the spectrum only at higher wavenumbers. The spectral form we observe suggests that such 
fine scales are prone to cluster into 2.5 Mm concentrations. 

Along similar lines, the imprint of structuring at still larger scales, such as by the 
network, manifests itself at lower wavenumbers, k < 0.2 rad/Mm. We find that actual quiet 
Sun spectra differ the most over these low wavenumbers, suggesting more pronounced day-to- 
day variation at these larger scales. Furthermore, any particular spectrum departs most from 
its exponential form over these low wavenumbers. This can be explained using fundamental 
elements distributed non-uniformly across the surface — concentrated, for example along 
network boundaries. This spectral range dictates the distribution of null points at great 
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heights, where very few ever occur. Thus the variation has virtually no effect on the total 
column we report. 

We performed our analysis on 562 different patches (300" x 300") of quiet Sun found at 
disk center over two different solar minima (1996-1998 and 2006-2007). The magnetograms 
were chosen to be free of active regions, plage or dominance of a single polarity. No other 
data was used to determine if the region was, for example, part of a coronal hole. Given their 
equatorial locations, however, it seems unlikely that much coronal hole data was included in 
our survey. We cannot therefore say whether the null point density applies to coronal holes 
as well as to equatorial quiet Sun. 

These findings place in a clearer context observations and models with putative rela- 
tionships to coronal null points. If the quiet Sun corona is being heated primarily through 
dissipation at null points, then the majority of this heat is supplied very low in the corona. 
The region between z = 1.5 Mm and z = 4.5 Mm contains three times as many nulls as the 
whole space above it; it would therefore probably receive a proportional share of the heating. 
It remains to be investigated whether the actual thermal structure of the quiet Sun corona 
is consistent with this distribution of heating, and is therefore consistent with heating by 
energy dissipation at null points. 

Along similar lines, if coronal jets were hypothesized to form anywhere a coronal null 
point occurs, then our results can be used in an ab initio prediction of their frequency. More- 
over, the distribution of their cusp heights should be somehow influenced by the distribution 
of null heights: short jets should out-number tall ones. Actual jets undoubtedly involve 
many factors, in addition to magnetic geometry, which influence the distribution of their 
properties. Still, it would be difficult to invoke null points as a requirement for jets if jets of 
a given height outnumbered null points at that level. A survey of coronal jets comprehensive 
enough to determine their frequency would thereby cast some light on the importance of null 
points in their construction. 



A. Estimating the modulation transfer function (MTF) 



To estimate the MTF, consider an ideal imaging system characterized by the same point 
spread function at every point on the image plane (i.e. an isoplanatic system). In this case the 
power spectral density (psd), S(k), of the image is the product of the actual psd, S(k ), and 
the square of the modulation transfer function (MTF), M(k) (IBorn and Wohill980l ). The 
MTF is the magnitude of the Fourier transform of the normalized point spread function. 
Generally, M(k), is a decreasing function of wave-number, causing the observed psd, S(k), 
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to be steeper than the actual one, S(k). 

If the same field is observed by MDI in both its high and low resolution modes the 
underlying psd, S(k), will be the same for both images. In this case 

SUk) jgg 

Shi(k) M^k) • 1 ) 

The MTF of the low resolution images will be dominated by the defocussing intended 
to reduce aliasing fr om spatial frequencies above the Nyquist limit k n = 1.6 rad/arcsec 



( iScherrer et a/.lll995l ). It is this contribution we can ascertain by cross comparison. 



One-minute high resolution magnetograms were obtained for 11 different intervals on 
five different dates throughout 2007 while MDI was in focus position 5. A single psd, S^k), 
was constructed for each time interval. In cases where several (typically 4) consecutive one- 
minute magnetograms were available from a given interval, the individual psd's were averaged 
together to produce S^k) with lower noise. The nearest 96-minute full disk low-resolution 
magnetogram was then used to compute S\ Q (k). The high and low resolution magnetograms 
were separated by six to 30 minutes (typically 10 minutes) strengthening the assumption 
that S(k) is the same for both observations. 

All psd's were computed as described in section 13.11 using spectral annuli of width 
Ak = 0.05 rad/arcsec, and the white noise floor was subtracted from each psd. The ratio 
of Equation flAlj) was computed over the range k < 1.5 rad/arcsec where both spectra were 
available and fairly free of noise. Full disk magnetograms have been re-calibrated by the 
MDI instrument team while high-resolution magnetograms have not. To compensate for the 
different calibration we divided the ratio by the mean of its lowest 8 spectral bins (k < 0.4 
rad/arcsec). 

A median MTF ratio is constructed by taking the median value of S\ Q (k) / Su(k) from 
all 11 time intervals for each spectral bin (see Figure [TTa) . It is fairly flat over the long- 
wavelength range and decays to low values by k > 1.2 rad/arcsec. Although the low reso- 
lution magnetograms extend out to k c — 1.8 rad/arcsec, it appears there is little sensitivity 
over that last third of that range. 

The estimated MTF ratio from Figure [TTa was fitted, over the range 0.1 rad/arcsec < 
k < 1 rad/arcsec, by a variety of empirical functions. Dashed and dotted curves show two of 
the more successful functions 



Mi/Mi = ^ , (A2) 



KIK = A'e-»' , (A3) 
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k [rod/orc.ec] 



Fig. 17. — The MTF and point spread function used for MDI. (a) The median of MTF 
estimates from 11 quiet Sun observations. The dark solid curve is the estimate from the 
ratio of Equation (lAlj) . Dashed and dotted curves are empirical fits. The triangle is the 
pre-launch value of the ratio quoted by Scherrer et al. (1995). (b) The point spread functions 
for the high resolution (dashed) and low resolution (solid) based on the MTF calculation. 
Sample pixel locations for each resolution are shown below the curves with crosses. 



respectively. The best fits were for parameter values £q = 1.04" and r = 0.88". The former 
curve better fits the flat inner region and the steeper fall-off beyond k ~ 0.6 rad/arcsec. There 
is evidence that the actual MTF falls even faster, perhaps even reaching zero, for k > 1.4 
rad/arcsec. Dividing the observed spectrum by expression (1A2I) is therefore a conservative 
correction at the very highest wave numbers. 

While this cross comparison provides the ratio of the two MTFs, it cannot give any 
information about either one alone. The most conservative correction to either spectrum is 
to assume the high -resolution mode is diffraction l imited and in perfect focus. In this case 
its MTF would be (IGhatak and Thyagarajan!ll978l ) 



M u (k) = M d {k) 



2 

7T 



COS 




k 2 



where k = 2nd/ = 
optical wavelength Am 



5.63 rad/ arcsec for a primary mirror of diameter, d 
= 6768 A JScherrer et a/ill995h . 



(A4) 



12.5 cm and 



Scherrer et al. (1995) reported that the intentional defocussing reduced the low reso- 
lution MTF at k — 0.8k n = 1.3 rad/arcsec to M — 0.17 from its diffraction-limited value 
of M d = 0.47. That is to say Mi /M hi = 0.36, under our assumption that M hi = M d . The 
square of this value is plotted as a triangle in Figure [T71 Our data point is consistent with 
this value, but the empirical fits are greater at this end of the spectrum. Using the empirical 
fits will then under-correct the highest wave-numbers. 
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